@ Rapid Calculation of Spacecraft Trajectories Using Efficient 
T aylor Series I ntegration 

Software greatly accelerates the calculation of spacecraft trajectories. 

John H . Glenn Research Center, Cle/eland, Ohio 


A variable-order, variable-step Taylor 
series integration algorithm was imple- 
mented in NASA Glenn'sSNAP (Space- 
craft N-body Analysis Program) code. 
SNAP is a high-fidelity trajectory prop- 
agation program that can propagate 
the trajectory of a spacecraft about vir- 
tually any body in the solar system. The 
Taylor series algorithm's very high 
order accuracy and excellent stability 
properties lead to large reductions in 
computer time relative to the code's 
existing 8th order Runge-Kutta 
scheme. Head-to-head comparison on 
near-Earth, lunar, Mars, and Europa 
missions showed that Taylor series inte- 
gration is 15.8 times faster than Runge- 
Kutta on average, and is more accurate. 
These speedups were obtained for cal- 
culations involving central body, other 
body, thrust, and drag forces. Similar 
speedups have been obtained for calcu- 
lations that include J2 spherical har- 
monic for central body gravitation. 
The algorithm includes a step size se- 
lection method that directly calculates 


the step size and never requires a re- 
peat step. 

High-order Taylor series integration al- 
gorithms have been shown to provide 
major reductions in computer time over 
conventional integration methods in nu- 
merous scientific applications. The objec- 
tive here was to directly implementTaylor 
series integration in an existing trajectory 
analysis code and demonstrate that large 
reductions in computer time (order of 
magnitude) could be achieved while si- 
multaneously maintaining high accuracy. 

This software greatly accelerates the 
calculation of spacecraft trajectories. At 
each time level, the spacecraft position, 
velocity, and mass are expanded in a 
high-order Taylor series whose coeffi- 
cients are obtained through efficient 
differentiation arithmetic. This makes 
it possible to take very large time steps 
at minimal cost, resulting in large sav- 
ings in computer time. The Taylor se- 
ries algorithm is implemented prima- 
rily through three subroutines: (1) a 
driver routine that automatically intro- 


ducesauxiliary variablesand setsup ini- 
tial conditionsand integrates; (2) a rou- 
tine that calculates system reduced de- 
rivatives using recurrence relations for 
quotients and products; and (3) a rou- 
tine that determines the step size and 
sums the series. The order of accuracy 
used in a trajectory calculation is arbi- 
trary and can be set by the user. The al- 
gorithm directly calculates the motion 
of other planetary bodies and does not 
require ephemeris files (except to start 
the calculation). The code also runs 
with Taylor series and Runge-Kutta 
used interchangeably for different 
phases of a mission. 

T his work was done by James R . Scott and 
Michael C. Martini of Glenn Research Cen- 
ter. Further information is contained in a 
TSP (seepagel). 

Inquiries concerning rights for thecommer- 
cial use of this invention should beaddressed 
to NASA Glenn Research Center, Innovative 
Partnerships Office, Attn: Steve Fedor, Mail 
Stop 4-8, 21000 Brookpark Road, Cleve- 
land, Ohio 44135. Refer to LEW-18445-1. 


@ Efficient Kriging Algorithms 

Goddard SpaceFlight Center, Greenbdt, Maryland 


More efficient versions of an interpo- 
lation method, called kriging, have been 
introduced in order to reduce its tradi- 
tionally high computational cost. Writ- 
ten in C++, these approaches were tested 
on both synthetic and real data. 

Kriging is a best unbiased linear esti- 
mator and suitable for interpolation of 
scattered data points. Kriging has long 
been used in the geostatistic and mining 


communities, but is now being re- 
searched for use in the image fusion of 
remotely sensed data. Thisallowsa com- 
bination of data from various locations 
to be used to fill in any missing data 
from any single location. 

To arrive at the faster algorithms, 
sparse SYMMLQ iterative solver, co- 
variance tapering, Fast Multipole 
Methods (FMM), and nearest neigh- 


bor searching techniques were used. 
These implementations were used 
when the coefficient matrix in the lin- 
ear system is symmetric, but not neces- 
sarily positive-definite. 

This work was done by Nargess 
M emarsadeghi of Goddard Space Flight Cen- 
ter. For further information, contact the God- 
dard Innovative Partnerships Officeat (301) 
286-5810. GSC-15555-1 


@ Predicting Spacecraft T rajectories by the WeavE ncke M ethod 

Lyndon B. Johnson SpaceCenter, Houston, Texas 


A combination of methods is proposed 
of predicting spacecraft trajectories that 
possibly include multiple maneuvers 
and/or perturbing accelerations, with 
greater speed, accuracy, and repeatability 
than were heretofore achievable. The 


combination is denoted the WeavEncke 
method because it is based on unpub- 
lished studies byjonathan Weaver of the 
orbit-prediction formulation of the noted 
astronomer Johann Franz Encke. Weaver 
evaluated a number of alternatives that 


arise within that formulation, arriving at 
an orbit-predicting algorithm optimized 
for complex trajectory operations. 

In the WeavEncke method, Encke's 
method of prediction of perturbed or- 
bits is enhanced by application of mod- 
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ern numerical methods. Among these 
methods are efficient Kepler's-equation 
time-of-flight solutions and self-starting 
numerical integration with time as the 
independent variable. Self-starting nu- 
merical integration satisfies the require- 


ments for accuracy, reproducibility and 
efficiency (and, hence, speed). Self-start- 
ing numerical integration also supports 
fully analytic regulation of integration 
step sizes, thereby further increasing 
speed while maintaining accuracy. 


This work was done by Jonathan K. 
Weaver of Johnson Space Center and Daniel 
R. Adamo of U nited Space Alliance For fur- 
ther information, contact theJSC Innovation 
Partnerships Office at (281) 483-3809. 
M SC -23802-1 


©An Augmentation of G -Guidance Algorithms 

This augmented algorithm can be used in small-body proximity operations utilizing model 
predictive control with a need for safety from surface-constraint uncertainty. 

N ASA's Jet Propulsion Laboratory, Pasadena, California 


The original G-Guidance algorithm 
provided an autonomous guidance and 
control policy for small-body proximity 
operations that took into account uncer- 
tainty and dynamics disturbances. How- 
ever, there was a lack of robustness in re- 
gards to object proximity while in 
autonomous mode. The modified G- 
Guidance algorithm was augmented 
with a second operational mode that al- 
lows switching into a safety hover mode. 
This will cause a spacecraft to hover in 
place until a mission-planning algorithm 
can compute a safe new trajectory. No 
state or control constraints are violated. 
When a new, feasible state trajectory is 
calculated, the spacecraft will return to 
standard mode and maneuver toward 
the target. The main goal of this aug- 
mentation is to protect the spacecraft in 
the event that a landing surface or obsta- 
cle is closer or further than anticipated. 
The algorithm can be used for the miti- 


gation of any unexpected trajectory or 
state changes that occur during standard 
mode operations. 

In order to have the G-Guidance algo- 
rithm detect an unsafe condition, it re- 
quired some modification. This modifi- 
cation provides a policy to safely 
maneuver the spacecraft between its cur- 
rent state and a desired target state while 
ensuring satisfaction of thruster and tra- 
jectory constraints, along with safety con- 
straints. In standard mode, this modifica- 
tion brings the spacecraft from its 
current position closer to its target state. 

I n safety mode, the algorithm maintains 
the spacecraft's current state at zero ve- 
locity. Since the safety mode is designed 
to be temporary, the destination location 
in this mode is also temporary, and once 
a new destination location is provided, 
the spacecraft returnsto standard mode. 

The G-Guidance algorithm uses both 
a planned trajectory (feedforward) and 


a control policy (feedback), along with 
sensors to monitor actual spacecraft 
state. The feedback is designed to en- 
sure that the spacecraft stays within a 
specified proximity to the feedforward. 
The feedforward is designed to achieve 
the goals of each mode: hover for safety 
mode and maneuver toward target for 
standard mode. By giving the spacecraft 
the ability to re-compute its trajectory 
on-the-fly in response to local condi- 
tions, minimization of fuel usage is pro- 
vided. The original G-Guidance algo- 
rithm provides robustness to uncertainty 
affecting the dynamics. The safety aug- 
mentation provides a form of state-con- 
straint robustness, which further miti- 
gates risk. 

This work was done by John M . Carson 
III and Behcet Acikmese of Caltech for 
NASA's Jet Propulsion Laboratory. For more 
information, contact iaoffice@jpl.nasa.gov. 
NPO-46452 


©Comparison of Aircraft Icing Growth Assessment Software 

The goal is to provide software that can predict ice growth under any condition for 
any aircraft surface. 

John H . Glenn Research Center, Cle/eland, Ohio 


A research project is underway to pro- 
duce computer software that can accu- 
rately predict ice growth under any me- 
teorological conditions for any aircraft 
surface. An extensive comparison of the 
results in a quantifiable manner against 
the database of ice shapes that have 
been generated in the NASA Glenn 
Icing Research Tunnel (IRT) has been 
performed, including additional data 
taken to extend the database in the 
Super-cooled Large Drop (SLD) regime. 
The project shows the differences in ice 
shape between LEWICE 3.2.2, Glen- 
nlCE, and experimental data. 


The Icing Branch at NASA Glenn has 
produced several computer codes over 
the last 20 years for performing icing 
simulation. While some of these tools 
have been collaborative projects, most 
have been developed primarily by one 
person, with some assistance by others. 
The state of computing has also 
changed dramatically in that time pe- 
riod. As these codes have grown in com- 
plexity and have been accepted by users 
as production icing tools, there has 
arisen a need for the developers to ad- 
here to standard software practices used 
to develop commercial software. 


The project addresses the validation 
of the software against a recent set of ice- 
shape data in the SLD regime. This vali- 
dation effort mirrors a similar effort un- 
dertaken for previous validations of 
LEWICE. Those reports quantified the 
ice accretion prediction capabilities of 
the LEWICE software. Several ice geom- 
etry features were proposed for compar- 
ing ice shapes in a quantitative manner. 
The resulting analysis showed that 
LEWICE compared well to the available 
experimental data. 

The effects of super-cooled large 
droplets in icing have been researched 
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